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Abstract 

The DiRC, a novel type of Cherenkov ring imaging device, is the primary hadronic particle identification 
system for the BaBar detector at the asymmetric B-factory Pep-II at SLAG. It is based on total internal 
reflection and uses long, rectangular bars made from synthetic fused silica as Cherenkov radiators and light 
guides. BaBar began taking data with colliding beams in late spring 1999. This paper describes the challenges 
for the DiRC in a higher luminosity environment and shows solutions to these challenges. 

1 Introduction 

The study of CP-violation using hadronic final states of the BB meson system requires the ability to tag the 
flavor of one of the B mesons via the cascade decay b ^ c ^ s, while fully reconstructing the final state of the 
other over a large region of solid angle and momentum. The momenta of the kaons used for flavor tagging extend 
up to about 2 GeV/c, with most of them below 1 GeV/c. On the other hand, pions from the rare two-body 
decays B'-* tt^tt^ (K^tt^) must be well-separated from kaons. They have momenta between 1.7 GeV/c and 
4.2 GeV/c with a strong momentum-polar angle correlation between the tracks (higher momenta occur at the 
more forward angles because of the cm. system boost) The Particle Identification (PID) system in BaBar[3] 
is located inside the calorimeter volume. Therefore, it should be thin and uniform in terms of radiation lengths 
(to minimize degradation of the calorimeter energy resolution) and small in the radial dimension to reduce the 
volume, hence, the cost of the calorimeter. 

The PID system being used in BaBar is a new kind of ring- imaging Cherenkov detector called the DiRC^ 
(the acronym DiRC stands for Detection of Internally Reflected Cherenkov light). It is designed to be able to 
provide excellent ir/K separation for all tracks from i?-meson decays from the pion Cherenkov threshold up to 
4.2 GeV/c. PID below 700 MeV/c exploits also the dE/dx measurements in the silicon vertex tracker and drift 
chamber. 



2 DiRC Concept 

The DiRC is based on the principle that the absolute values of angles are maintained upon reflection from a flat 
surface. Figure shows a schematic of the DiRC geometry that illustrates the principles of light production, 
transport, and imaging. The radiator material of the DiRC is synthetic, fused silica in the form of long, thin bars 
with rectangular cross section. These bars serve both as radiators and as light pipes for the portion of the light 
trapped in the radiator by total internal reflection. Fused, synthetic silica (Spectrosil[S]) is chosen because of its 
resistance to ionizing radiation, its long attenuation length, large index of refraction, small chromatic dispersion 
within the wavelength acceptance of the DiRC, and because it allows an excellent surface polish on the bars[H]. 
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Figure 1: Schematic of the DiRC fused silica radiator bar and imaging region. 



In the following, the variable 9c is used to designate the Cherenkov angle, 4>c denotes the azimuthal angle 
of a Cherenkov photon around the track direction, and n represents the mean index of refraction of fused silica 
(n = 1.473) within the wavelength acceptance of the DiRC (300 nm to 600 nm). The Cherenkov angle is given 
by the familiar relation cos^c = ^/{nP) with f3 — v/c, the velocity of the particle v and the velocity of light c. 

For particles with [3^1, some photons will always lie within the total internal reflection limit, and will 
be transported to either one or both ends of the bar, depending on the particle incident angle. To avoid 
instrumenting both ends of the bar with photon detectors, a mirror is placed at the forward end, perpendicular 
to the bar axis, to reflect incident photons to the backward, instrumented end (see Figure [T]l . 

Once photons arrive at the instrumented end, most of them emerge into a water-filled expansion region, 
called the standoff box. A fused silica wedge at the exit of the bar reflects photons at large angles relative to 
the bar axis. It thereby reduces the size of the required detection surface and recovers those photons that would 
otherwise be lost due to internal reflection at the fused silica/water interface. The photons are detected by an 
array of densely packed photomultiplier tubes (PMTs)|3IH|, each surrounded by reflecting light catcher cones [^l 
to capture light which would otherwise miss the active area of the PMT. The PMTs are placed at a distance of 
about 1.2 m from the end of the bars. The expected Cherenkov light pattern at this surface is essentially a conic 
section, where the cone opening-angle is the Cherenkov production angle modified by refraction at the exit from 
the fused silica window. 

The DiRC bars are arranged in a 12-sided polygonal barrel (see Figure EJ. Because of the beam energy 
asymmetry (at the T{4S), Pep-II collides 9 GeV electrons on 3.1 CeV positrons), particles are produced 
preferentially forward in the laboratory. To minimize interference with other detector systems in the forward 
region, the DiRC photon detector is placed at the backward end. The bars are placed into 12 hermetically sealed 
containers, called bar boxes, made of aluminum-hexcel panels. Dry nitrogen gas flows through each box, and is 
monitored for humidity to ensure that the bar box to water interface remains sealed. Each bar box contains 12 
bars, for a total of 144 bars. Within a bar box the 12 bars are optically isolated by a ~ 150 /Ltm air gap between 
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Figure 2: Exploded view of the DiRC mechanical support structure. The iron magnetic shield is not shown. 

neighboring bars, enforced by custom shims made from aluminum foil. 

The bars are 17 mm thick, 35 mm wide, and 4.9 m long. Each bar is assembled from four 1.225 m pieces 
that are glued end-to-end; this length is the longest high-quality bar currently obtainable [Bl El- The bars are 
supported at 600 mm intervals by small nylon buttons for optical isolation from the bar box. Each bar has a 
fused silica wedge glued to it at the readout end. The wedge, made of the same material as the bar, is 91 mm 
long with very nearly the same width (33 mm) as the bars and a trapezoidal profile (27 mm high at bar end, and 
79 mm high at the light exit end). The bottom of the wedge (see Figure^ has a slight 6 mrad) upward slope 
to minimize the displacement of the downward reflected image due to the finite bar thickness. The 12 wedges in 
a bar box are glued to a common 10 mm thick fused silica window which provides the interface and seal to the 
purified water in the standoff box. 

The standoff box (see Figure is made of stainless steel, consisting of a cone, cylinder, and 12 sectors of 
PMTs. It contains about 6000Z of purified water. Water is used to fill this region because it is inexpensive and 
has an average index of refraction (n « 1.346) reasonably close to that of fused silica, thus minimizing the total 
internal reflection at silica-water interface. Furthermore, its chromaticity index is a close match to that of fused 
silica, effectively eliminating dispersion at the silica-water interface. The iron gusset supports the standoff box. 
An iron shield, supplemented by a bucking coil, surrounds the standoff box to reduce the field in the PMT region 
to below 1 G[n. 

The PMTs at the rear of the standoff box lie on an approximately toroidal surface. The distance from 
the end of the bar to the PMTs is ~ 1.17 m. Each of the 12 PMT sectors contains 896 PMTs (ETL model 
9125[7||H1) with 29 mm diameter, in a closely packed array inside the water volume. A double o-ring water seal 
is made between the PMTs and the vessel wall. The PMTs are installed from the inside of the standoff box 
and connected via a feed-through to a base mounted outside. The hexagonal light catcher cone is mounted in 
front of the photocathode of each PMT which results in an effective active surface area light collection fraction 
of about 90%. 

The DiRC occupies 80 mm of radial space in the central detector volume including supports and construction 
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Figure 3: Display of an e+e^ M^A*^ event reconstructed in BaBar. with two different time cuts. On the left, 
all DiRC PMTs with signals within the ±300 ns trigger window are shown. On the right, only those PMTs with 
signals within 8 ns of the expected Cherenkov photon arrival time are displayed. 

tolerances with a total of about 19% radiation length thickness at normal incidence. The radiator bars cover a 
solid angle corresponding to about 94% of the azimuth and 83% of the cm. polar angle cosine. 

3 DiRC Performance 

The performance of the DiRC is influenced by two factors. First, information from the PMTs has to be separated 
into background and signal. Second, the 9c resolution of the signal has to allow for the separation of different 
particle species. 

Figure 121 shows a typical di-muon event (e+e^ fi'^i-i^). In addition to the signals caused by Cherenkov 
light from the muon tracks, about 500 background signals can be seen in the 600 ns readout window centered 
around the trigger. This background is dominated by low energy photons from the Pep-II machine hitting the 
standoff box. 

The time-to-digital converter (TDC) chip 12 used in the DiRC data readout is designed such that a dead 
time of about 5% occurs at an input rate of 250 kHz. Some care in machine tuning is required to stay under a 
limit of 250 kHz per tube. To monitor this rate, one PMT in each sector is read out via a scaler. Figure 0] shows 
the maximum scaler rate as a function of the Pep-II luminosity during data taking in 2000 and 2001. In March 
2000 the accelerator operated at a peak luminosity of lO'^'^ cm~^ s~^. Figure0Ja) shows, at a value corresponding 
to only one third of the design luminosity, the PMT rates reached a level that caused noticeable dead times. 
Due to those findings, lead shielding was installed in the summer of 2000 around the beam line components near 
the backward end-cap. Initially, localized shielding was added in the form of lead bricks which were stacked 
around the beam pipe and in front of a large quadrupole. This shielding significantly improved the background 
situation so that noticeable TDC dead times were reached only at 2.5 x 10'^^ cm~^ s~^ (Figure EJb)). During 
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Figure 4: Maximum scaler rates measured for three different accelerator and shielding configurations in March 
2000 (a), October 2000 (b) and August 2001 (c). 



the shutdown in January 2001, the localized lead brick shielding was replaced by an engineered, homogeneous 
lead shielding of 5 cm to 8 cm thickness covering the inside radius of the standoff box and is easily removable 
to facilitate access to the central detector and beam line components. As is shown in figure QJc), the maximum 
scaler rates at luminosities of 4.2 x 10^'^/ cm~^ are well below the level that would cause TDC inefficiencies. 

During the shutdown following the 2001-2002 run, we plan to replace the TDCs with a faster version with 
deeper buffering that is designed to have a deadtime of less than 5% at 2.5 MHz input rate, suitable for 
luminosities of at least 10'^^ cm~^ s~^. 

The DiRC is intrinsically a three-dimensional imaging device. Photons are focused onto the phototube 
detection surface via a "pinhole" defined by the exit aperture of the bar, so that the photon propagation angles 
ax and ay can be measured in two-dimensional space, where x (bar width) and y (bar thickness) are the directions 
transverse to the bar axis. The travel time of the photon down the bar is also related to the photon propagation 
angle (az) with respect to the bar axis. Imaging in the DiRC occurs in all three of these dimensions, by recording 
the time of the PMT signal. As the track position and angles are known from the tracking system, these three 
a angles can be used to determine the two angles {9c,4'c)- This over-constraint on the angles is particularly 
useful in dealing with ambiguities and high background rates. 

The single photon resolution can be calculated from the geometrical, the chromatic and the transport term. 
The geometrical uncertainty is, for the pinhole optics of the DiRC, given by the standoff box size, the bar 
and the PMT size. With 7.2 mrad, it is the single biggest contribution to the single photon resolution. The 
chromatic uncertainty originates from the fact that the photons are produced with different wavelengths and, 
therefore, have different cone opening angles 9c- Using the range in photon wavelength accepted by the PMTs, 
the chromatic effect contributes to the resolution with 5.4 mrad. The final uncertainty comes from imperfections 
of the bars and the glue joints. This transport term is approximately 2 mrad to 3 mrad. 

Combining the three contributions, a single photon resolution of uec-i — 9-5 mrad is expected for the DiRC. 
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Figure 5: The difference between the measured and expected Cherenkov angle for single photons, ^9c,-y for 
single muons in fJ,~^fi~ events. The curve shows the result of a fit of two Gaussians to the data. The width of 
the narrow Gaussian is 9.6 mrad. 

Figure |S1 shows the single photon angular resolution obtained from di-muon events. There is a broad background 
of less than 10% relative height under the peak that originates mostly from track-associated sources, such as 6 
rays, reflections off the glue-fused silica boundaries, and combinatorial background 13 . The width of the peak 
translates to a resolution of about 9.6 mrad, in good agreement with the expected value. 

The resolution on the track Cherenkov angle is defined by a correlated term and the single photon resolution 
which scales with the number of photoelectrons: 



where Np^ is the number of photons detected per track. 

The number of photoelectrons varies between 16 for small values of cos 6'track at the center of the barrel and 60 
at large values of cos^track as is shown in Figure El This variation is well reproduced by Monte Carlo simulation 
and can be understood from the geometry of the DiRC. The number of Cherenkov photons varies with the path 
length of the track in the radiator, it is smallest at perpendicular incidence at the center and increases towards 
the ends of the bars. In addition, the fraction of photons trapped by total internal reflection rises with larger 
values of cos^track. This increase in the number of photons for forward going tracks is a good match to the 
increase in performance required at larger momentum. 

With the present alignment, the track Cherenkov angle resolution for di-muon events is shown in Figured 
The width of the fitted Gaussian distribution is 2.4 mrad compared to the design goal of 2.2 mrad. 

The measured time resolution is 1.7 ns, close to the intrinsic 1.5 ns transit time spread of the PMTs. This 
resolution is used to efficiently distinguish background from signal photons but is not sufficient to improve the 
Cherenkov angle resolution. 

Figure IHl shows an example of the use of the DiRC for particle identification. The Ktt invariant mass spectra 
are shown with and without the use of the DiRC for kaon identification. The peak corresponds to the decay of 
the particle. 

The FID performance of the DiRC has been studied with a sample of pions and kaons, selected kinematically 
using — > K^TT^ decays from inclusive D* production. 
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Figure 6: Number of detected photons versus track polar angle for reconstructed tracks in di-muon events 
compared to Monte Carlo simulation. The mean number of photons in the simulation has been tuned to match 
the data. 




Figure 7: Resolution of the reconstructed Cherenkov polar angle per track for di-muon events. The curve shows 
the result of a Gaussian fit with a resolution of 2.4 mrad. 
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Figure 8: Invariant Ktt inclusive mass spectrum with and without the use of the DiRC for kaon identification. 
The mass peak corresponds to the decay of the D'^ particle. 



The tt/K separation power of the DiRC was defined as the difference of the expected Cherenkov angles for 
pions and kaons, divided by the measured track Cherenkov angle resolution. As shown in FigurelHl the separation 
between kaons and pions at 3 GeV/c is about 4.4cr, within 10% of the design goal. 

The efficiency for correctly identifying a charged kaon that traverses a radiator bar and the probability to 
wrongly identify a pion as a kaon are also determined from the inclusive D* sample and are shown as a function 
of the track momentum in Figure 1101 for a particular choice of particle selection criteria. The kaon selection 
efficiency and pion misidentification, integrated over the K and tt momentum spectra of the D* control sample, 
are 97.97 ±0.07% (stat. only) and 1.83 ±0.06% (stat. only), respectively. 

A report on the operational experience with the DiRC detector can be found in [TJ . 



4 DiRC in Higher Luminosity Environment 

In a higher luminosity environment, the accelerator-induced background will cause the main challenge. From a 
set of special runs taken in February 2002, the background rate as a function of the beam currents /her and 
Iler in the high and low energy ring, respectively, and the luminosity C can be parameterized as 

kHz kHz kHz ^ 

R - 13^/HER + 18^ W + 10 ^^33 (2) 

For the conditions at the beginning of this year, the current terms dominate but for higher luminosities, the 
corresponding term will contribute significantly. The anticipated limit of the luminosity Pep-II will be able to 
provide without major upgrade is 4 x 10'^^ cm~^ s~^. With some upgrades to the shielding and the TDC upgrade. 
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Figure 9: DiRC n/K separation versus track momentuni measured m.D^ ^ K 7r+ decays selected kinematically 
from inclusive D* production. 
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Figure 10: Efficiency and misidentification probability for the selection of charged kaons as a function of track 



momentum, for a particular choice of particle selection criteria, 
kinematically from inclusive D* production. 
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Figure 11: New design of DiRC standoff region using focusing optics. 



the DiRC will be able to reliably work in this environment. At higher luminosities, such as expected after major 
Pep-II updates or for SuperBABAi?,, more changes will be necessary. 

A solution to high background rates is to reduce the size of the standoff box. However, without modifications 
to the design, this would lead to an increase in the geometric resolution term and would cause a considerable 
decrease in the DiRC performance. 

A modified design, as shown in figure ^2 employs smaller sized photodetectors as well as focusing optics 
effectively removing the bar size uncertainty in the focusing plane. These two enhancements can be balanced so 
that the size of the standoff region can be reduced while keeping a similar or improved geometrical resolution. 

In addition, smaller photodetectors with a better transit time spread will improve the timing resolution of 
the detector. This can be used to further reduce the background by using tighter time cuts and to reduce the 
chromatic uncertainty. The travel time of a photon is the difference between the time the particle hits the bar 
and the arrival time of the photon. It depends on the path length and the group velocity in the transit medium. 
The group velocity itself depends on the refractive index. Thus, the travel time and path length yield a measure 
of the wavelength of the photon. Figure El shows a calculation demonstrating that a timing resolution in the 
order of 100 ps is needed for a significant improvement in the chromatic uncertainty. 

Reducing the range of wavelength accepted in the detection system will reduce the range of Cherenkov cone 
opening angles contributing and, thus, reduce the chromatic uncertainty. The drawback of this solution is the 
loss in the number of photons detected per track causing a reduced suppression of the single photon resolution 
term. These two effects can be balanced in order to reach the optimal total Oq resolution. 

With a timing resolution in the order of 100 ps, the ratio of the track to detector distance over the path 
length will provide an additional measure of the photon angle. In some regions of phase space, especially at 
small dip angles, its resolution will be competitive to the angle resolution from the position of the PMT hit. 

A more detailed discussion can be found in |ll5j . 

One of the photon detector candidates is the flat panel PMT H-8500 by Hamamatsu'16. . It has 8x8 square 
pads on a size of 50 mm x 50 mm. This leads to a high area efficiency and packaging density. The timing 
resolution is expected to be in the range of 100 ps to 200 ps. The disadvantage of the currently available model 
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Figure 12: Achievable wavelength resolution as a function of photon wavelength assuming 100 ps and 500 ps 
timing resolution. 

is its low gain in the order of 1.6 x 10^. 

Studies of the timing resolution on an early pre-production version were encouraging. The time was measured 
at two thresholds for each signal in order to perform a time walk correction. The corrected time distribution 17 is 
shown in figure lT^ A fit to a double Gaussian plus a second order polynomial describes the data well. The narrow 
Gaussian has a width of 125 ps. This value is not yet corrected for the light pulse width and instrumentational 
effects. In addition, the double threshold method limits the data to PMT signals with a pulse height above the 
higher threshold. This removes some fraction of low signal events. 

Currently, studies of the efficiency uniformity of the PMT area are in process and results will be available in 
near future. 

5 Improvements for Future DiRC Detectors 

The design described in the previous section will be sufficient to provide the required resolution of the DiRC for 
possible upgrade scenarios. However, for the SuperBABAR environment with a design luminosity of 10'^^ cm~^ s~^, 
the projected occupancy in the tracking system makes building a completely new detector necessary. 

This allows to reconsider the radiator and light guide design. Synthetic fused silica is still the material of 
choice as it fulfills all requirements in terms of transparency, radiation hardness, achievable uniformity in material 
and optical finish. In addition, it has a small radiation length which is benifical for a calorimeter outside of the 
DiRC system. 

Increasing the thickness of the radiators will lead to an increased number of photons per track but it will also 
increase the amount of material in front of the calorimeter, so the thickness of the radiator can not be changed 
significantly. 

Increasing the width of the radiators, such as using one plate instead of several bars, has advantages in 
the geometrical term of the single photon resolution. In a plate geometry, there are less side bounces and the 
path length for a photon measured in a given PMT with a known track position changes drastically with each 
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Figure 13: Time resolution of the Hamamatsu H-8500 PMT. The data was recorded using two different thresholds 
to allow for time walk correction. No other correction has been applied. 



additional side bounce. Therefore, the hit time can be used to determine the number of side bounces and because 
of this the exit position of the photon at the instrumented bar end can be inferred. This is shown in figure [HI 
where the full line and the dashed line have significantly different path lengths. 

There are a few examples in the following on the expected resolution for the different setups. A detailed 
discussion of the geometries can be found in JS]. Please note that the resolutions are given in terms of the 
photon exit angles a at the end of the bar or plate. These are different to the Cherenkov angle resolution oq^ . 
However, there is a strong correlation between (Jq,^ and a^c while the correlation between CTq^ and crgc weaker. 

Current BaBar Dirc 

The current DiRC setup has round PMTs with 29 mm diameter and rectangular bars of size 35 mmx 17.5 mm. 




Figure 14: Plate geometry. 
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The resolution on the photon exit angle at the edge of the bar is 



1 1 , 1 



+ (pixel) ^ YTto^V 12^^^ 16^^^ " ^^'^ 

and 

(T„^ « 7.6 mrad (4) 

The interface fused sihca — water magnifies the angles thus improving the resolution by approximately 10% 
leading to cr^^^^ = 9.5 mrad/6.9 mrad. 

Focusing optics for BaBar like DiRC 

A focusing optic in y with a pixel size of 6 mm x 6 mm and a corresponding standoff distance of i = 250 mm 
leads to 

» V<^x (pixel) » 2b{)mm V 12 " ^'^ 

similar to the current ay resolution. The resolution in x will get worse as focusing in both dimensions would 
complicate ambiguity resolution. This can be compensated by the reduced chromatic uncertainty. 

Focusing optics with plates for SuperBABAR 

In order to improve the resolution, a smaller pad size of 2 mm in y direction and 6 mm in x direction can be 
used. The resolution can be described by 



l^/^IP^ « ^ V 1^(2 mm)^ = 2.3 mrad (6) 



in y direction and in x 



~ \ (track) + a2 (pixel) « ^ (4 mm)2 + 1 (6 mm)2 = 4.3 mrad (7) 

In the last equation L is the length of the region where the photon can expand without side bounces. This length 
varies depending on the number and positions of the side bounces. 

Together with the improved chromatic resolution, this leads to a significant improvement in the single photon 
resolution compared to the current BaBar Dirc. 

The major uncertainty not directly related to the DiRC design is the correlated uncertainty on the position 
of the track at the entry into the radiator. In order to improve the particle separation, it is not sufficient to 
improve the intrinsic DiRC uncertainty only. A possible solution is to add a tracking device outside the DiRC 
volume. This would also allow for finding particle conversions and decays inside the radiator material. 

A further improvement is to enlarge the angular coverage of the particle identification system by using an end- 
cap device. The Dirc covers the barrel region only. Additional challenges for an end-cap device are the limited 
amount of space in the end-cap, the higher radiation levels close to the beam-pipe, and that the readout has to 
be done inside the magnetic field which makes standard PMTs unusable. Interesting alternative photodctcctors 
have been suggested but they still need to be studied in detail before an evaluation on the feasibility of an 
end-cap Dirc can be performed. 
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6 Conclusions 



The DiRC is a novel particle identification system used for the first time in the BaBar environment. It is 
performing close to the design and has a significant influence on most physics analysis. 

Increasing the luminosity beyond the BaBar. design values will result in challenges for the DiRC system mainly 
due to increased background levels of low energy photons from Pep-II. Improvements in shielding and in the 
maximal data taking rate will allow for running the current detector up to the anticipated Pep-II limits. Beyond 
that, reducing the size of the standoff box will be necessary. In order to keep the Chcrenkov angle resolution 
on the design level, the standoff box needs to be replaced by focusing optics with smaller sized photodetectors. 
Timing resolution in the order of 100 ps is vital to reduce the chromatic uncertainty and can effectively separate 
signal from background photons. Research on candidate devices has started and first tests have shown promising 
results. 

At luminosities in the order of 10*^^ cm^^ s^^ as envisioned for SuperBABAR, a completely new DiRC can be 
built. Fused silica plates instead of bars and even smaller PMT pad sizes along with an improved tracking will 
allow for significant improvements in the Cherenkov angle resolution. 

Enlarging the angular coverage of a DiRC particle identification system will require an end-cap device. The 
challenges for this device include photodetection inside of the magnetic field and the limited amount of space 
prohibiting large expansion volumes. 
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